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Abstract: In this article we propose a novel test for statistical anisotropy of the CMB 
Ar(n= (6^,0)). The test is based on the fact, that the Galactic foregrounds have a 
remarcably strong symmetry with respect to their antipodal points SI iv ^ — n, n = 
(0,0) and 5^2 := n ^ n, n = TT — 0, with respect to the Galactic plane, while the 
cosmological signal should not be symmetric or asymmetric under these transitions. 
We have applied the test for the octupole component of the WMAP ILC 7 map, by looking 
at aa^i and a3^3, and their ratio to a3^2 both for real and imaginary values. We find 
abnormal symmetry of the octupole component at the level of 0.58%, compared to Monte 
Carlo simulations. By using the analysis of the phases of the octupole we found remarkably 
strong cross-correlations between the phases of the kinematic dipole and the ILC 7 octupole, 
in full agreement with previous results. 

We further test the multipole range 2 < / < 100, by investigating the ratio between the 
/ + m = even and l-^m — odd parts of power spectra. We compare the results to simulations 
of a Gaussian random sky, and find significant departure from the statistically isotropic 
and homogeneous case, for a very broad range of multipoles. We found that for the most 
prominent peaks of our estimator, the phases of the corresponding harmonics are coherent 
with phases of the octupole. We believe, our test would be very useful for detections of 
various types of residuals of the foreground and systematic effects at a very broad range 
of multipoles 2 < / < 1500 - 3000 for the forthcoming PLANCK CMB map, before any 
conclusions about primordial non-Gaussianity and statistical anisotropy of the CMB. 
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1. Introduction 

Since the release of the COBE temperature map of the Cosmic Microwave Background 
(CMB), and after the Wilkinson Microwave Anisotropy Probe (WMAP) ([1], [2]), the mod- 
ern cosmology has been dramatically changed. Two fundamental hypotheses, the existence 
of the Dark Matter and the Dark Energy in nature, have found experimental confirmation 
in the CMB data. Moreover, the power spectrum of the temperature anisotropy and po- 
larization of the CMB is in excellent agreement with theoretical predictions of the simplest 
models of inflation, providing a remarkably important connection between the physics of 
the modern Universe and the very beginning of the cosmological creation of the matter. 
However, right after the COBE mission, some anomalies of the CMB sky attracted very 
serious attention, starting with the low amplitude of the quadrupole [3], the quadrupole- 
octupole alignment with planarity and lack of correlations at > 60^ [4], [5], [6], the phase 
correlations between different multipoles [7], [8], [9], [10], [11], the existence of the cold 
spot [12], [13], [14], the dominance of the power of odd multipoles over even one [15], [16], 
[17], [18], [19], [20] and [21], the ecliptic and Galactic north-south asymmetry of the power 
spectrum, [22], etc (see for review [23]. 

In this article, we present a new test for statistical asymmetry and possible non- 
gaussianity of the CMB, based on an investigation of the symmetries of the CMB sky with 
respect to different directions on the sphere. This symmetry test will reflect the properties 
of the Galactic foregrounds and possible systematic effects. The idea of the test is based on 
the assumption that the primordial CMB signal, unlike any other non-cosmological signal 
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(like the galactic synchrotron, free- free, dust emission and galactic and extragalactic point- 
like sources), should have no particular symmetries or asymmetries of angular distribution 
on the sky. This is what one would expect for a completely random, chaotic distribution. At 
the same time, for instance, the galactic foregrounds clearly have a well defined symmetry 
of the temperature distribution with respect to their antipodes, and slightly broken, but 
still "visible" symmetry with respect to the Galactic plane [24]. The instrumental noise, 
the residuals of the calibrations and possible other effects of systematics have a preferable 
symmetry in the direction of the north and south ecliptic poles (in Galactic coordinates). 

All these symmetries manifest themselves in the coefficients ai^j^ of the spherical har- 
monic decomposition of the signal AT(0, 0) trough the even and odd multipoles (the sym- 
metry SI = h ^ — fi, fi = {0^(j)) with respect to the antipodes), trough the symmetry 
5^2 = n ^ n with respect to the Galactic plane, where n = (tt — 0,0), and trough the 
symmetry 53 = fi ^ m, where m = {0,27v — (/)). 

Note, that the octupole component is the most powerful tail of the CMB power spectrum. 
Before our analysis, the peculiarity of the octupole was widely discussed in connection 
with quadrupole- octupole alignment [3], [4], [6]. In this article, we take a closer look at the 
various ratios between the a^^^-values making up the octupole, unconnected with any in- 
vestigation of the quadropole, and find a significant deviation from a random distribution. 
Further, we create power spectra from the real and imaginary part of the a^^^-parameter, 
where the sum is over even and odd values of / + m only {D^{1) and D~{1) respectively). 
Then we test the ratio between real and imaginary Z)^(/)-values, for a large range of mul- 
tipoles. The new element of our analysis, which was never discussed before, is that the 
symmetry of / = 7 mode of the ILC 7 map is peculiar at the level of 3 events from 10^ 
realizations for / + m = odd. For l^m — even the most impressive result is connected with 
/ = 14, also with 3 events from 10^ realizations, while for / = 38 we have the corresponding 
probability about 0.6%. 

It would be worth to note, that our analysis is based on the ILC 7 map, which is 
contaminated by the point-like sources (Galactic and extra-galactic, residuals of the diffuse 
foregrounds, uncertainties of the antenna beam and possible effects of systematics). This 
is why the ILC map can not be used for evaluation of the CMB power spectrum without 
implementation of the mask (for instance, the KQ75 WMAP mask). After that, the anal- 
ysis of the a/^^-coefficients can not be done, due to a very strong coupling between the 
coefficients, induced by the mask. However, our method is especially useful for estimating 
the degree of contamination of the ILC map - in combination with standard methods for 
determining the power spectrum C{1) from the masked sky - in order to check out possible 
sources of peculiarities of C{1). We believe, that this method would be especially useful for 
the ongoing PLANCK mission, where the ILC map would be applicable to a very broad 
range of multipoles, compared to the WMAP range. An important point is, that for very 
high /, actually, we do not have any theoretical predictions about potentially dangerous 
zones, which need to be masked. 

There is one more important implementation of our method, related to the non- 
Gaussianity test (the so called fni- approach), allowing one to constrain different models 
of inflation. Quadratic corrections to the linear theory of perturbations from inflation (lo- 
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cal fni models), are characterized by a coupling between low and high multipoles, which 
provides fingerprints of this particular type of non-Gaussianity in the bi-spectrum, and in 
higher order moments (see for review [25]). To assess the f^i approach with maximal pre- 
cision, it is clear that we have to detect all possible sources of non-Gaussian contamination 
of the CMB map, amongst which the point sources seems to be a major component. 

The outline of the paper is the following. In section 2 we introduce the various symme- 
tries of the CMB sky, and analyze the ratio between various values of ai^^n foi" the octupole 
(Section 3). In addition, in Section 3 we will show that the symmetry test is closely 
connected with the phases of the coefficients of the spherical harmonic decomposition, in- 
dicating the most peculiar components. Further, in Section 4 we introduce the pathfinder 
of peculiar multipoles, based on the symmetry test, created from even and odd values of 
/ + m, and finally we apply that method to the ILC 7 map. In section 5, we summarize 
the results of simulations and comparison with the data, given by our estimators. 

2. The symmetry test 

The temperature fluctuations on the CMB sky can be decomposed into spherical harmonics 
in the following standard way: 

Ar(n) ^Y.Y. ^hmYlAi^) = -7= J2J^^ai,m=o)Pi{cose) + 

1=0 m=-l 1=2 ^ 

+ ^ E E / (^^ + iy-/")' p,'"(cosg) X [SR(a;,^)cos(m^) - S(a;,^) sin(m0)] 

(2.1) 

where ai^^n is the coefficient of decomposition, and n = (0, (/)), with and (/) being the polar 
and azimuthal angles on the sky respectively. 3? and denote the real and imaginary parts 
of the a^^^-coefficients, P^^(cos0) are the associated Legendre polynomials, and Pi{cos0) = 
Pl^{cos0). For the CMB signals AT(n) on the sphere, one can define the symmetric 
AT^(n) and anti-symmetric AT~(n) components: 

AT(n) = AT+(n) + AT"(n), (2.2) 

where all inversions have been given with respect to the origin of the polar system of 
coordinates, which is a standard basis for estimation of the a/^^-coefficients according to 
Eq(2.1), and 

I m=—l 

Ar(n)-Ar(-n) „ 
Ar-(n) = ^ = E E ai,mP-{l)YUn), 

I m=—l 

(2.3) 
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Figure 1: Top left. Classification of the symmetries of the CMB with respect to the antipodal 
points (SI), the Galactic plane (S2) and (j) ^ 27r — (j) (S3). Top right. The map of a signal with all 
multipoles and m-modes. Then, from the second from the top row (left) and down to the bottom 
right panel we show the maps for / = even^ I = odd^ l-\-m = even, l-\-m = odd, l-\-m = even, = 0, 
l-\-m = even, 5Re = 0, / + m = odd, = 0, and l-\-m = odd, 3?e m = 0. All in Galactic coordinates 
and Mollweide projection. 

and P+(/) = cos2(f ), P-{1) = sin2(f ). 

By definition, AT+(n) and AT~(n) are orthogonal in the sense that the averaged 
product over the whole sphere vanishes (i.e. (Ar+(n)AT~(n)) = 0). The power spectrum 
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of AT(n) is given by 

^(^) = ^ E = ^^(0 + C-(0, (2.4) 

m 

where = C{l)P^{l) and = C{l)P-{l). In other words, and are 

associated with the power spectrum of even and odd multipoles respectively. 

The area near the Galactic plane contain a very powerful concentration of the Galactic 
diffuse foregrounds (synchrotron, free- free and dust emission), and is highly contaminated 
by the Galactic point-like sources. The Galactic plane is perpendicular to the z-axis and 
centered at z = in the Galactic coordinate system and, as we mentioned above, is the basis 
of a type of symmetry. The coordinate inversion from n = (0, 0) to n = (tt — 0, (j)) give us a 
symmetric ATs(n) = lSTs{^) oi" an anti-symmetric ATa(n) = — Ara(n) signal with respect 
to that plane. Since Y/^^(n) = (—1)^+^ Y/^^(n), one can see that the symmetric signal AT^ 
corresponds to / + m = even and AT^ corresponds to / + m = odd. Taking into account 
Ar+ = Ar+ + AT+ and AT" = AT^ + AT^", we can see that the most symmetric 
part of the signal AT (with respect to the origin and the plane at z = 0) corresponds 
to AT+ with / = even^ I ^ m — even, and the most anti-symmetric part is AT" with 
/ = odd and I + m — odd. We have summarized all the mentioned symmetries with the 
illustration in Fig.l. Note that all a/^^ coefficients were found from the polar system of 
coordinates, centered at the Galactic center. Then, by selecting different multipoles (for 
instance, only even /, or only odd), we plot the corresponding maps in Galactic coordinates 
and in Mollweide projection. 

3. The ILC 7 octupole. General properties. 

As an illustration of the different kind of symmetries in the CMB sky, we plot the octupole 
component of the ILC7 map in Fig. 2-3. Note that the octupole is the most powerful anti- 
symmetric component with respect to the inversion fi — > — fi on the sky. For the octupole 
(i.e. / = 3), even and odd m number corresponds to odd and even values of / + m. From Eq. 
2.1, one can see that symmetry (/) ^ 27r — 0, = const for even and odd m is associated with 
real and imaginary part of the spherical harmonic coefficients (i.e. 3?[a;^^] and ^[ai^jyi]). 




Figure 2: The map of the octupole in Galactic (left) and ecliptic (right) coordinates. 
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Figure 3: Components of the ILC7 octupole: On the left and the right side, we show the octupole 
components, where only aim of even and odd l-\-m are retained. At the same time, we have retained 
both real and imaginary part of aim (top), real part of aim (middle) and imaginary part of aim 
(bottom). 



Coming back to Eq(2.1, we can see that different symmetries or asymmetries of the 
octupole component are connected with corresponding components of as^rn coefficients, 
which now play a role of weighting coefficients for different symmetric or asymmetric terms 
(see table 1). For a statistically homogeneous and isotropic random Gaussian field the 
distribution of the phases of as^m foi" m = 0, 1, 2, 3 is uniform in the interval [0, 27r], while 
the amplitudes |a/^^| are distributed according to the Rayleigh probability density function. 

3.1 Symmetry estimators as an example of Cauchy distribution. 

We now investigate the a/^^'s of the octupole. It is clear, that the value of a/^^ that 
best conforms to the criteria for asymmetry is as^2^ where / = odd and I ^ m — odd. 
Naturally very many values of a/^^ meet this criteria, but a3^2 is the one with most power, 
as mentioned above. We now wish to test the ratio between the a3^2-niultipole, and its more 
symmetric neighbors aa^i and a3^3 (where / + m = even), for both the real and imaginary 
value. To do this, we introduce the following estimators: 



^(Q^3,i) 
^^(^3,2)' 



as = 



^(0^3,3) 
^^(^3,2) 



3^(a3,i) 
5ft(a3,2)' 



/33 = 



3^(^3,3) 
3«(a3,2) 



(3.1) 
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If all the real and imaginary parts of 
the estimators follow a random Gaussian 
process, the parameters «i,«3,/3i and /33 
follow a Cauchy distribution function. This 
particular type of random process is char- 
acterized by the probability density func- 
tion: 



/ = 3, 


m 


3?e(a3,m) 




m — 





-6.479e - 03 





m — 


1 


-1.219e-02 


2.0265e - 03 


m — 


2 


2.199e-02 


5.907e - 04 


m = 


3 


-l.lTle - 02 


3.355e - 02 



A 



{x — xo)^ + 7^ 



(3.2) 



Table 1: 5Re and parts (in mK) of the mul- 
tipole coefficients as^m for the ILC 7. The cor- 
responding power for / + m = even is Ceven = 
2.0210~^mK'^ ^ and for / + m = odd we have 
Codd = 7.5110-^mi^2 



where A is the normalization constant, xq 
is the location parameter, and 7 is the probable error. The Cauchy distribution is one 
of the examples of a random distribution without mean value, variance or higher order 
moments. However the probability P(x > X) to get some corresponding values of x > X 
is still defined, and it is given by the integral in Eq(3.3). 

For the random Gaussian field, when nominators and denominators in Eq(3.1) are normal- 
ized to ^/0.5C{l), we have A^I/tt, xo = and 7=1, and thus the result in Eq(3.3) can 
be simplified significantly. 



P{x > X) = A 



X 



dx 

(x - x^y + 7^ 



A 

7 



TT 



— tan 



X - Xq 

7 



ttX' 



X > 1 (3.3) 
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Figure 4: Value of the estimators from eq. 3.1, compared to a histogram of 10.000 Monte Carlo 
simulations. The blue line corresponds to ai, the red line is a^. The green line show Pi the vertical 
black line (somewhat obscured by the green) corresponds to Ps- Note that the y-axis is in log. The 
smoothed black line corresponds to the Cauchy distribution. 
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We have performed the test on the WMAP 
ILC 7 year data, and compared with 10.000 Monte 
Carlo simulations of the CMB sky. Note that we 
run numerical simulations in order to assess the 
level of fluctuations for finite sample and X <^ 1. 
One more remark is connected with the prop- 
erties of the Cauchy distribution, given by esti- 

mators from Eq(3.1). Since ah the estimators ^^^^^ ^, ^^^^^ ^ ^ ^ 

are based on the ratios of (potentially) Gaus- \a^l\asl\pil\P3\ for the WMAP octupole 
sian variables, the big values of them seems to and the 10.000 Monte Carlo simulations 
be very unlikely, as well as small values X ^ 0. (MC) in comparison with P{x > X) from 
The point is that inversions \ai\ and Eq(3.3)(Th). 

lAI ~^ transform small values of the pa- 

rameters to big one, leaving the Cauchy distribution in form of Eq(3.2) unchangeable. The 
results are plotted in fig. 4, and in table 2. 

We see a big deviation for the 0:3 parameter, with only 58 out of 10.000 simulations 
having a larger value, corresponding to 0.58%, while from Eq(3.3) we have corresponding 
probability P{x > X = 56.8) = 0.0056. For the ai parameter, we have 9.4%, and for 
the two f3 parameters, the deviation is not significant. Thus it seems, that the imaginary 
part of the as^2 multipole carries little weight, especially compared to the a^^s multipole. 
It is clear, that the octupole exhibit a symmetric behavior, and contributions from the 
foreground seems to be one of the most likely explanations. However, in the next section we 
will show that this problem could be more complicated, than the very natural assumption, 
that only the foreground component is responsible for the peculiar symmetry of the WMAP 
octupole. 

3.2 The WMAP 7 octupole. The phases. 

As it was pointed out in Section 2, the S3 symmetry reflect directly the properties of real 
and imaginary parts of the ai^^n coefficients. In application to the octupole component, the 
corresponding estimators ai,as,Pi and were designed to evaluate different symmetries 
or asymmetries of the WMAP octupole. However, for a given value of m the ratio between 
imaginary and real parts of the a/^^ is nothing, but the phase of the corresponding /,m 
harmonic, defined as 

- tan-' f^^^\ (3.4) 

According to [7], [8], the phases of the ai^ra coefficients reflect directly the morphology of 
the AT-map, and for the octupole component they are connected with q^i,q^3,/3i and /^a 
as follows: 

tan(^3,i) ^ tan(^3,3) ^ 

Q^i — ri — ^3 — -r, — rPs (3.5) 
tan('^/;3^2) tan(^3^2) 

The analysis of the phases is very useful for investigating possible cross-correlations 
between two or more signals, a^^^ and //^^, where fi^^a corresponds to the foreground com- 
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ponent (residuals of the Galactic diffuse foreground and extra-galactic point-like sources), 
for instance, the WMAP V or W bands. At the same time, the multipole vectors approach 
[4], [6] clearly detect the coupling between the WMAP quadrupole and octupole and kine- 
matic (non-cosmological ) dipole. This result was never tested before by analysis of phases 
of the octupole component, as presented below. In Fig. 5 we show the phases of the WMAP 
ILC7 octupole and NILC5 octupole from [26] versus the phases of the kinematic dipole. 




-1.0 I , , , , \ , , , , \ , , , , \ , , , , I 

-1.0 -0.5 0.0 0.5 1.0 

cos«-) 



Figure 5: The phases for NILC5 (black) and ILC7 (blue) octupole. The solid line is for ips,i = 2.976 
rad, the dash line is for iIjs,2 = 0.0268 rad. and the dash-dotted line is for i/j^^s = 1.907 rad. The 
red solid line is for dipole q = 0, and the red dash line is for = ^^/^ phases. 

Due to periodicity of the phases within the interval — tt, tt, the estimators of cross- 
correlations between different signals are based on the trigonometric moments ^ = 
cos('0/^^ — ^^;/^^/)(see [8] for review). Thus for the WMAP 7 octupole and the phases of 
the kinematic dipole, the corresponding estimators are : 

Glj = -0.986, G3;2 = 0.99964, G]'^^ = 0.9440. 

One can see that the phase correlations of the ILC 7 octupole and the kinematic dipole 
is especially strong for 3, 2 component, practically at the same level of significance, as 
for 0:3 test. Does it mean, that the major source of peculiar symmetry of the octupole 
is connected to the residuals of the dipole substructure? Answering that question, we 
would like to stress, that in Fig. 5 the phases of the NILC5 octupole are slightly different 
from the corresponding phases of the ILC 7 octupole. The NILC approach differs from 
the WMAP ILC method by different evaluation of the CMB-foreground coupling. Taking 
under consideration the templates derived by the WMAP team for the synchrotron, free- 
free, and dust emission, we have synthesized the sum of all these components for the V 
band, in order to obtain the phases $3,rn of the octupole component, listed below: 

$3 = 0, $3,1 = -0.0342, $3,2 = 1.5126, $3,3 = -2.47629 
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and further: cos('03^i - $3^1) = -0.991, cos('03^2 - ^3,2) = 0.0849 and sin('03^2 - ^3,2) = 
0.9963, and cos('^3^3 — $3,3) = —0.323. Thus, from the analysis of the ipi^m and ^i^rn cross- 
correlations, we may conclude that 3, 1 component of the octupole anti-correlate with 3, 1 
phase of the foreground, and 3, 2 component of the octupole is orthogonal to the 3, 2 
component of the foreground. This is why, we believe, the simplest assumption that only 
the foreground or only the kinematic dipole are responsible for the abnormal symmetry of 
the ILC 7 octupole, does not reproduce the whole picture of the possible contamination of 
the primordial signal. 

Lets illustrate this idea by an analysis of the following model of contamination. Our 
assumption is, that the WMAP octupole component contains the primordial signal cs^rm 
contaminated by the residuals of the foregrounds fs^m and residuals of the kinematic dipole 
as follows: 



1 

tt3,m cs^m + ^ ^(3, m|l, m')di^m' + M/3,m, m = 1, 2; (3.6) 

where: L(3,m|l,m^) is the linear shift operator, which takes the dipole components to 
3, 2 and 3, 3 components of the octupole, and fi is the "ILC-residuals of the foreground" 
coupling parameter. For the imaginary part of m = 2 component of the octupole from 
Eq(3.6) we get: 

1^3,21 sin(^3,2) = |c3,2|sin(^3,2) + I ^ L(3, m|l, m')di,m^ J + m|/3,2| sin($3,2), 

(3.7) 

Since sin('03^2) — 0, from Eq(3.7) one can get: 



|c3,2| sin(^3,2) - -"^rn I ^ L(3, m|l, m')di^rn' j - m|/3,2| sin($3,2) (3.8) 

Thus, in framework of linear model of the ILC contamination, the primordial component of 
the octupole is fully determined by the residuals of the foregrounds and kinematic dipole. 
This result clearly illustrate the importance of investigation of the symmetries for different 
components of the ILC signal, if we are interested in a more accurate estimation of the 
statistical properties of the CMB signal, especially for high multipoles /. 



4. S2 -pathfinder of peculiar multipoles. 

In the previous section we have shown that the implementation of the symmetry test is 
very informative in application to the analysis of the anomalies of the low multipole tail 
of the CMB power spectrum. However, for high multipoles / » 10, the number of possible 
permutations of real and imaginary parts of the corresponding a/^^ coefficients grows very 
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rapidly. This is why for this range of multipoles it would be essential to propose the 
quick search test, based on the above mentioned symmetries SI — S2. As the basis of this 
test, we compute the power spectrum, using the real part and the imaginary part of ai^ 
respectively, where the sum only includes either even or odd values of / + m. 



C{1) = Z)(0+e + D{l)-Re + mL + mjm^ 

I 



1 



^« 21 + 1 



m=l 



(4.1) 



m=l 



where G^{l^m) — cos^( ^^^^^^ ) and G~{l^m) — sin^( ^^^^^^ ). We are going to estimate 
the ratio — D{1)^^/ D{1)^^ for / + m = even components and = ^i^ReZ-^^^^Jm 
l^m — odd in order to investigate the statistical isotropy and Gaussianity of CMB signal. 
Effectively, we are thus only comparing power spectra composed of either symmetric (+) 
or asymmetric (— ) terms, with respect to the S2 symmetry. 

We have tested the 7-statistic for the WMAP 7 ILC data, until / = 100, and compared 
them with 1000 Monte Carlo simulations. Below, in Fig. 6, we plot the results of our 
analysis for both even and odd / + m components of C{1). 

As is evident, we have several peaks for the 7^^, which is most noticeable at / = 14, 
where only 3 events out of 1000 realizations have values as high as the WMAP data. For 
that multipole the phase of 14,4 component is '^14,4 = —3.061 c:^ — tt, and '^14, 10 = 2.920, 
which indicate the most planar components. 

Other prominent peaks for are 
found at / = 38 and / = 97, where only 
6 and 3 events have value as large as 
WMAP data respectively. For the 7^", 
we see a similar tendency: at / = 7 we 
have 3 events higher than that of the 
WMAP data. Another notable peak is 
at / = 47, where we have 4 events higher 
than the WMAP data. 
The 7^~ have some very prominent peaks 
for / < 10 (at / = 3 for instance, the 
ratio is in the order of 10^, and for / = 7, 
the ratio is around 10^). These peaks 

are to high to include in the plot above, as they would have dwarfed the other values 
completely, and as we wanted the ratios to be easily comparable, we decided against a 
logarithmic y-axis. We summarize in Fig. 7-8 the images of the corresponding maps for the 
given values of / and all m, for the most prominent peaks from Fig. 6. One can clearly see 
the corresponding symmetries or asymmetries of these maps with respect to the galactic 
plane. For illustration of the morphology of these maps, in Fig. 9 we plot the phases of 
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Table 3: 3?e and parts (in mK) of the multipole 
coefficients aj rn 

for the ILC 7. 
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Figure 6: Left panel. The ratio (7^"^) between C{l)Re and C{l)im for even / + m: The red hne 
shows the result from the ILC7 map, the black horizontal line shows the unit value, and the black 
dots correspond to the distribution of 1000 Monte Carlo simulations. Right panel. The same as 
left, but for jf. 




Figure 7: The images of the WMAP ILC 7 components with / = 14 (left), / = 38 (middle), and 
/ = 97 (right panel) 



the ai^m for / = 7 and / — 14, combining I ^ m — even^ and I + m — odd modes. Even 
without special analysis, one can see that phases of / + m = even modes are coherent with 
the phases of the octupole. 
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Figure 8: The images of the WMAP ILC 7 components with / = 7 (left), / = 47 (middle), and 
/ = 68 (right panel) 




-i.ol , , , , \ , , , , \ , , , , \ , , , , I -i.ol , , , , \ , , , , \ , , , , \ , , , , I 

-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

cos((') cos((') 



Figure 9: The phases of the WMAP ILC 7 components with / = 7 (left), and / = 14 (right panel). 
The blue line indicate m = component, the black lines correspond to / + m = odd^ the red dash 
lines are for / + m = even. The green lines show the phases of the octupole. 

The peculiar behavior of the multipoles, detected by the 7+ and 7"- tests, allow us to 
look closely at the corresponding symmetry estimators for each. In table 3 we show the 
corresponding real and imaginary parts of the a/^^-coefficients for / = 7 peak. From this 
table, it is clearly seen that for / = 7 the ratio a^^i = 11.08. This parameter is 

an analogue of the 0^3- parameter for the octupole and it corresponds to the S^S-symmetry. 
Taking under consideration that the corresponding probability for that parameter is given 
by Eq 3.3 with X = 11.08, we get P{a7^i) 2^ 0.028. The analogue of the /3-parameter of 
the octupole is ^7^4 = - 21.69, with P(/37,4) ^ 0.015. Thus, for / = 7 the most 

established type of symmetry is S2, with an additional S3 component. 

5. Conclusion 

In this paper, we have introduced 3 symmetries on the sphere of the sky: one with respect 
to the antipodal points (SI), one with respect to the galactic plane (52), and one where 
we rotate in the ^-direction {S3). The a/^^-coefficients, would either be symmetric or 
asymmetric with respect to the three symmetries, depending on the value of /, / + m and 
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the real and imaginary components of the a/^^-components. Therefore we have introduced 
a symmetry test, based on the ratio between the symmetric and asymmetric ai^m- 

We have tested it for the WMAP7 octupole in particular, as it is the most powerful 
of the asymmetric multipoles, and the results were odd at the level of 0.58% for 0^3, and 
9.4% for the ai compared to 10000 Monte Carlo simulations. The quadru- and octupole 
have previously been extensively investigated in connection with the problem of alignment, 
but this result is only based on the octupole, unconnected with the quadropole. We have 
shown, that the octupole carry anomalously little power in the imaginary part of the most 
asymmetric component, as^2^ especially compared with the imaginary part of the more 
symmetric a^^s- 

We have further introduced a symmetry test, based on the ratio between power spectra 
of real and imaginary parts of the a^^^, created from sums over even and odd values of 
/ + m only (7^"^ and respectively). We tested this for the WMAP7 data, in the range 
2 < / < 100, and compared with 1000 Monte Carlo simulation. We found notable deviations 
from the simulations at / = 14 (0.3%), / = 38 (0.6%) and / = 97 (0.3%) for 7+, and at 
/ = 7 (0.4%) and / = 47 (0.4%) for 7". 

We have tested all mentioned anomalies by the phase analysis and find that all these 
symmetries and asymmetries are closely related to the correlations of the phases. We 
have investigated the coupling of the octupole phases with phases of dipole and foreground 
and confirm remarkable correlations of the octupole with the kinematic dipole. At the 
same time we note, that residuals of the foregrounds could play a significant role in the 
formation of the peculiar symmetry of the octupole. We believe, that our symmetry test 
would be very useful for estimating the quality of the separation of the primordial CMB 
from non-cosmological signals, making the analysis of the primordial non-Gaussianity more 
sensitive. 
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